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Structure formation models with a cosmological constant are successful in explaining large-scale 
structure data, but are threatened by the magnitude-redshift relation for Type la supernovae. This 
has led to discussion of models where the cosmological 'constant' decays with time, which might 
anyway be better motivated in a particle physics context. The simplest such models are based on 
scalar fields, and general covariance demands that a time-evolving scalar field also supports spatial 
perturbations. We consider the effect of such perturbations on the growth of adiabatic energy 
density perturbations in a cold dark matter component. We study two types of model, one based 
on an exponential potential for the scalar field and the other on a pseudo-Nambu Goldstone boson. 
For each potential, we study two different scenarios, one where the scalar field presently behaves 
as a decaying cosmological constant and one where it behaves as dust. The initial scalar field 
perturbations are fixed by the adiabatic condition, as expected from the infiationary cosmology, 
though in fact we show that the choice of initial condition is of little importance. Calculations 
are carried out in both the zero-shear (conformal newtonian) and uniform-curvature gauges. We 
find that both potentials allow models which can provide a successful alternative to cosmological 
constant models. 
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I. INTRODUCTION 



Ever since the demise of the Standard Cold Dark Matter model, which proved unable to simultaneously match 
the COBE observations of cosmic microwave background anisotropics, the abundance of rich galaxy clusters and the 
shape of the galaxy correlation function, it has been fashionable to study variations on this basic theme. A popular 
alternative is to lower the matter density, thus shifting the epoch of matter-radiation equality in a desirable direction. 
To keep the possibility that standard models of inflation, generating a spatially-flat Universe, might be responsible for 
the seed perturbations, the most popular version of this model is to introduce a cosmological constant A to maintain 
the spatial flatness. This has proven a compelling framework within which it is possible to understand the formation 
and evolution of large-scale structure in the Universe ||l|,|| . 

The favoured models, until recently, had the energy density J^a in the cosmological constant lying in the range 0.5 
to 0.7. As well as giving a good fit to large-scale structure observations, these models have support from two further 
sources. If one assumes the baryon density in the Universe is that predicted by standard big bang nucleosynthesis 
1^, and that large galaxy clusters provide a representative sample of the universal ratio between baryons and the 
total amount of matter in the Universe then the matter density must be significantly below one. Further, the 
need to have the age of the Universe, t^, exceeding the age of the globular clusters in our galaxy suggests that I^a 
should be as large as possible. Until recently the required value for was usually around 14 Gyrs but a revised 
Cepheid distance scale due to new distance estimates by the satellite Hipparcos has brought this down to something 
more like 12 Gyrs 1^,0. Given that most measurements (even taking into account Hipparcos's revision of the Cepheid 
distance scale) seem to suggest that the present value of the Hubble parameter, h, is at least 0.6 (in the usual units 
of 100 kms-^Mpc"^) H and perhaps even larger in a flat Universe we would then need J7a > 0.3 (0.55 if 
to > 14 Gyrs). 

Unfortunately, these models have been dealt a serious blow by the preliminary results from the Supernova Cosmology 
Project jl^], which attempts to determine the magnitude-redshift diagram of Type la supernovae. They place a 95 
per cent confidence upper limit of J^a < 0.5 for flat Universes. This is significantly stronger than earlier limits from 
the galaxy velocity distribution ||ll[], galaxy outflows from voids and the statistical analysis of the frequency of 
gravitational lensing of high-redshift quasars [|3|, all coming in around ilA < 0.7. While the ACDM model remains 
viable with these smaller JIa values, this is seen as much less attractive because, as in the critical-density case, the 
required matter density is well above that given by direct observation. 

A way out of this dilemma is to move to cosmological models where the cosmological constant is substituted by a 
dynamical quantity which decays with time ||l^|23|| . Within these models it is possible to relax the constraints resulting 
both from the frequency of gravitational lensing of high-redshift quasars and from Type la supernovae | 
While in some ways this is clearly a regressive step, introducing more freedom into the model, it can also be argued 
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that such a situation may be more natural on particle physics grounds. For example, there is the well-known difficulty 
within quantum field theory to understand the very small vacuum energy density, /Xvac = (0.003 eV)^f2A, required by 
a cosmological constant. If not strictly zero, due to some yet unknown cancellation mechanism, one would expect 
/ivac to be between 50 and 120 orders of magnitude larger A decaying cosmological constant term would be a 
simple way of reconciling a very large vacuum energy density early on in the Universe with an extremely small one 
at present. 

Some authors have simply assumed more or less ad hoc decay laws for the cosmological constant term [ p^ . By 
comparing predictions of the models with observations it was then hoped that the correct decay law could be recovered, 
which would then shed some light on the possible physical process behind the decaying cosmological constant term. 
However, a credible mechanism for obtaining such a term already exists, which is to assume the existence of a scalar 
field presently relaxing towards the minimum of its potential (see e.g. p5|-p3[). Scalar fields are not only predicted 
to exist by some particle physics theories that go beyond the Standard Model, but are also the most plausible engine 
behind a possible inflationary period in the very early Universe [ p7| . The overall dynamics of the Universe in the 
presence of a relaxing scalar field, and its consequences for several classical cosmological tests, has been studied in 
detail by various authors |24y2§|] . 

Our aim in this paper is to study the effect of spatial perturbations in the cosmological 'constant' term given by 
such a field, in particular on the growth of the matter perturbations. When one has the standard constant A term, 
there is no possibility of any perturbations in it. However, as soon as one permits any form of time variation, general 
covariance immediately implies that it must be able to support spatial perturbations. Despite this, presumably 
because of the development of this line of research from the original constant case, with few exceptions ]l9|~p3[| most 
authors have not looked at the possible effect of spatial perturbations in the background value of the scalar field on 
the growth of perturbations in the matter distribution. 



II. EQUATIONS AND INITIAL CONDITIONS 

We will assume that the background space-time contains an ideal fluid and a scalar field. The equation of state of 
the ideal fluid, relating its background pressure, Pj, to its background energy density, /j,^, is = (7 — 1)M7j where 7 
is a constant. The background energy density and pressure associated with a minimally coupled real scalar field with 
potential V{(j)) are given by 

W = ^'^' + V(0) ; p^ = ^^^-V{^). (1) 

The overdots represent derivatives with respect to coordinate time t. 

The evolution of these background quantities is described by the Friedmann equation. 
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together with the two energy conservation equations 



A7 = -iHjfi^ , (3) 
4> + 3H^ + = , (4) 

for the ideal fluid and the scalar field respectively. Here H = a/a is the Hubble parameter, a the cosmic scale factor, 
and ,0 represents a derivative with respect to </>. 

We shall carry out a fully relativistic treatment of the perturbations, using a formalism based on a series of papers 



by Hwang |29-33|. As always, the equations describing the evolution of the perturbations are very complicated and we 
have relegated their discussion to our two Appendices. A crucial question when studying perturbations is the choice 
of gauge. We shall consider two choices. Mainly we shall use the zero-shear gauge (ZSG) [also known as conformal 
newtonian or longitudinal], but we shall also solve the equations in the uniform-curvature gauge (UCG) as a check 
on the accuracy of the numerical calculations. 

The perturbation equations, described in the Appendices, relate the perturbed part of the metric variables, a 
(perturbed part of the lapse function), ip (perturbed part of the spatial curvature), x (perturbed part of the shear) 
and K (perturbed part of the expansion scalar), to the perturbed part of the matter variables, e = e-y + (perturbed 
part of the total energy density), m — zu-y + vj^ (perturbed part of the total pressure) and ^E" = ^'-y + '^^ (perturbed 
part of the total energy density flux, or total fluid four- velocity, depending on the frame chosen). In the ZSG x 
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is chosen to vanish, while in the UCG vanishes. In all, we will need to numerically integrate a s yste m o f seven 
simu ltan eous first-order ordinary differential equations, formed by the background e quati ons, E qs. ( Bl) to ( |B3| ), where 
Eq. (^^) g ives /i -y, an d the four perturbation equations resulting from either Eqs. ( Bl^ ) and ( Big| ) (in the ZSG), or 



Eqs. (P2q) and (|B21|) (in the UCG). 

We must now specify the initial conditions required in order to solve this system. We shall do that by providing 
the initial values for the quantities e-y, ^-y, Scf) and 54)' . The fi rst t wo giv e the initial values for and (ZSG) (x and 
X' in the UCG) by means of Eqs. (|b|) and (^ [Eqs. ( pT^ ) and (|B^) in the UCG]. 

We make two requirements on the initial conditions. The first is to choose the initial value of so that only the 
growing mode of the solution for the evolution of an ideal fluid in an Einstein-de Sitter Universe is present (see the 
exact analytical solutions in Table 1 of Ref. ^^), as this is the situation we would expect in the later phase of the 
evolution of the Universe if the density perturbations in the ideal fluid were produced in the very early Universe. 

Secondly, we will assume that the energy density perturbations in the ideal fluid and the scalar held are related by 
the adiabatic condition Q. Standard models of inflation Q always give this type of perturbation, as there is only 
one dynamical degree of freedom during inflation, and so it is by far the most natural choice to make. The entropy 
perturbation should therefore vanish. 



together with its flrst time derivative, 
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where 



and 



The quantities 
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represent the internal entropy of each component, 
conditions we obtain the initial values of 56 and 5 



In the case of the ideal fluid we have e-, = 0. From the adiabatic 



~ (37M7)[2^>0-fc'<^'/(3a)] 



(10) 
(11) 



We also need to specify the initial values of the background variables H, (j) and (/)' which determine the cosmological 
model. The initial value of fi-y is obtained by the requirement that the Universe has critical density. The other 
three degrees of freedom for the initial background conditions are fixed by requiring specific present values of the 
cosmological quantities h, to and Sl!^. Note that for some V{(l>), and given the required values for h and ri°, it may 



not be possible to obtain the desired value for t^. 
that can be reached for a given ¥{(/)). 



In general, for fixed h and fl^, there will be a maximum value 



III. THE CHOICE OF SCALAR FIELD POTENTIAL 



Before we go on to perform the numerical integration of the perturbation equations in both the ZSG and the UCG, 
we need to specify what type of ideal fluid we will consider and the shape of the scalar fleld potential. We are mainly 
interested in the growth of matter density perturbations during the matter dominated era, and there are good reasons 
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to believe that most matter in the Universe has neghgible intrinsic velocity, i.e. it is cold. In this paper we will assume 
that the ideal fluid has no pressure, that is 7 = 1. 

With relation to the type of scalar field we should consider the issue is not as clear. Many different scalar fields 
have been proposed, most of them with the specific aim of producing an inflationary expansion phase in the very early 
Universe, though some also originate from attempts at extending the Standard Model of particle physics. 

We will consider two different scalar field potentials, an exponential potential of the form 

1/(0) =Uoexp(-/30), (12) 

and the potential associated with a pseudo-Nambu-Goldstone-boson (PNGB) field js^, 

y(0) = m4 [cos(0//) + 1] . (13) 

Both these potentials have been extensively studied, within the context of power-law inflation |3^] and natural inflation 
psi respectively. The PNGB field has also been proposed explicitly as the most natural candidate for a presently- 
existing minimally- coupled scalar field [|l^,^ . 

In both cases we have two degrees of freedom, and in principle one should explore the full 2-dimcnsional parameter 
space defined by them. However, our main objective in this paper is to draw attention to the importance of taking 
into account the possibility of spatial perturbations in a cosmological scalar field when one is assumed to exist, due 
to their influence on the growth of matter density perturbations. Therefore, in this paper we will only consider two 
sets of values for the constants associated with each of the two potentials. 

In two of the models thus obtained, one for each potential, the scalar field presently behaves like a slowly-decaying 
cosmological constant. In these models, we will choose the initial values of the background variables H, cj) and 4>' , 
and the constants associated with the potentials, so that we end up with h — 0.6, il^ — 0.4, and an age for the 
Universe of to = 14 Gyrs. In the other two models, again one for each potential, the scalar field starts behaving like 
non-relativistic matter, scaling as a~^, at a redshift close to 100. Therefore, the age of the Universe in these models 
is pretty much the same as if the Universe was always matter dominated. In order to obtain an acceptable age we 
lower h slightly to h — 0.55 to give to — 12 Gyrs. 

Our numerical integration of the perturbation equations will begin at a redshift of z — 1100, roughly at electron- 
photon decoupling, and end at the present time. One reason for this choice is that through the COBE satellite 
measurement of the amplitude of cosmic microwave background anisotropics we have good knowledge of the amplitude 
of energy density perturbations existing at the horizon scale at this redshift ]39|] . The main reason though is that 
for reasonable parameter values we are well into matter domination and the effects of radiation on the matter power 
spectrum of density perturbations have already run their course. At this redshift, the power spectrum is well described 
by the cold dark matter transfer function, for example as parametrized by Bardeen et al. |^,^. We are thus able 
to analyze the effects of the scalar field on the matter power spectrum, without them being concealed within the full 
Boltzmann code machinery. We cannot however make predictions for the full microwave anisotropy power spectrum. 



A. The exponential potential 

For the exponential potential there is a particularly interesting situation, which we will call EXPl, where the relative 
energy densities of the scalar field and the ideal fluid remain constant with time, thus implying — (7 — l)/^^, after 
a transitional period. The energy density associated with the scalar field is then a fixed fraction, 247rG'7//3^, of the 
total energy density in the Universe [ [l6p^j2^ . This scaling solution of the cosmological background equations is one 
of two attractor points for the system, the other being the well-known power-law inflation solution. The former is 
the one which is attained if the potential is steep enough. As homogeneous perturbations around the scaling solution 
typically have complex eigenvalues, the system usually approaches the attractor point through oscillations in the 
relative energy densities of the scalar field and the ideal fluid. It has recently been studied by Ferreira and Joyce , 
though for different parameters than the model we will look at. 

That the Universe eventually reaches the scaling solution with Jim — 0.4 demands that (3 — vIOttG. We choose 
the scalar field energy density at the start of the simulation to be much smaller than that of the scaling solution; this 
is not necessary, though it does have to be true much earlier at nucleosynthesis [ ^2|j2^ . This restricts the possible 
combinations for the values oi 4> at z — 1100 and Vq. The value (j> takes at z = 1100 is a matter of definition and 
we set it to zero. We will assume cj) to be extremely small at z = 1100, thus implying that Vo = (0.025 eV)"*. With 
this choice of parameters the scalar field begins to contribute significantly to the total energy density in the Universe 
by a redshift of about 70. In accordance with the discussion above, this model therefore gives the same age as the 
Einstein-de Sitter case, and to make it large enough we choose h = 0.55. 
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The second situation we consider, denoted EXP2, is where the system is presently entering a period of power-law 
inflationary expansion, analogous to the exponential expansion which would arise if a cosmological constant is present 
in the Universe. To obtain our desired values (stated above) for Jl^, h and to implies roughly that (3 = \/28ttG. On 
the other hand, these constraints only very weakly define the initial splitting of the energy density associated with the 
scalar field aX z — 1100 into its potential and kinetic parts. We will again arbitrarily assume that (f){z = 1100) ~ 0, 
thus implying Vb = (0.0025 eV)"*. 



B. The PNGB potential 



In the case of the PNGB potential we will also examine two distinct possible evolution histories for the Universe. In 
the model we will cah PNGBl, where M = 0.04 /i^/^ eV and / = 1.25 x 10^^ GeV, the scalar field starts to roll down 
its potential at a redshift of about 100, and is presently coherently oscillating at the bottom of the potential, behaving 
dynamically as cold matter. This model, like EXPl, gives the same age as a Universe containing only matter, so again 
we take h to be 0.55, giving an age for the Universe of 12 Gyrs. The constraint il"^ = 0.4 implies that 0(z = 1100) = / 
and 4i{z = 1100) ~ 0. 

In model PNGB2, where M = 0.003 /i^/^ eV and / = l/%/8^ = 2.4 x 10^^ GeV, the scalar field begins moving 
down its potential at around z = 1, and has not yet reached the bottom of the potential by the present time. The 
values of M and / in this model coincide with those chosen for a more detailed analysis by Frieman et al. in their 
paper dealing with PNGB motivated dynamical cosmological constant models.]^ Imposing the boundary condition 
= 0.4 leads to (j){z = 1100) = 1.75/ and 0(z = 1100) ~ 0, while h = 0.6 and to = 14 Gyrs are obtained by simply 
choosing the correct value for H aX z = 1100. 



IV. RESULTS 



We will work in the ZSG, and display the fluid density perturbation 5^ = e;^//i^. We arbitrarily take its initial 
value at redshift 1100 to be 10~^; as the equations are linear, and our main results only show power spectra relative 
to one another, the initial magnitude of the density perturbation is irrelevant. Note that the density perturbations 
given in different gauges coincide on scales well within the horizon, but otherwise, though uniquely defined as long 
as the gauge is given, do not coincide. Gare must therefore be taken in interpreting any long-wavelength behaviour. 
Most of the literature uses the comoving gauge if the large-scale power spectrum is shown. 

We integrated the background and perturbation equations both in the ZSG and the UCG using the NAG Fortran 
Library Routine D02CJF, which is based on a variable-order, variable-step Adams method, for values of k in the 
range of 10~^ to 1 /iMpc^^. We did this not only for the four scalar field models chosen in the previous section, but 
also for two cases where no scalar field was assumed present, which will serve as comparison. In one the Universe has 
critical density, while in the other the Universe is flat with JIq = 0.4, thus implying the presence of a cosmological 
constant. 

In Fig. 1 we plot the present amplitude obtained for 5'!^{k) (for simplicity henceforth omitting the superscript and 
always meaning this gauge) for our four models, relative to that found in the case of a flat Universe with ilo = 0.4 
(with h in this comparison model adjusted to match that of each of the four scalar field models we consider). The 
main integration runs were done in the Z SG, b ut when the same runs were performed in the UCG, with €^{z = 1100) 
calculated from e^{z — 1100) using Eq. (B28) and then e^^{z = 0) from el^(z — 0) using Eq. ( B31 ), the difference in 
the present value of was less than one per cent over the whole range for k in the four models. 

Although we believe that Eqs. (jlffl) and (O) should be used to obtain the initial values for 5(j) and (50', in accordance 
with the adiabatic relations Eqs. (0) and (0) one expects from inflationary generated perturbations, we also looked 
at what occurs when arbitrary initial values for these quantities are assumed. After performing several integration 
runs for a variety of initial values for 5<t) and Scj)' , we reached the conclusion that the evolution of the energy density 
perturbations in the ideal fluid is almost (to a few per cent) independent of the initial values one considers for these 



^In a more recent paper. Coble et al. considered a slightly different model with M = 0.005 eV and / = 1.885 x 10^* GeV, 
where the scalar field starts moving down the PNGB potential at z ~ 10 and is presently already oscillating around its minimum. 
These values were chosen so that the model yields = 0.4 and h — 0.7. This implies an age for the Universe only slightly 
above 10 Gyr. We prefer the PNGB2 model over theirs. If one chose to obtain h = 0.55, and thus an age close to 12 Gyr, one 
would end up with fi^j = 0.5. 
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FIG. 1. This shows the present amplitude of fractional perturbations in the energy density of the ideal fluid in the ZSG 
for our four different models, relative to the one obtained in the case where no scalar field is present and the Universe is flat 
with fio = 0.4. The full and dashed lines represent respectively the two exponential potential and the two PNGB models. The 
arrow indicates the inverse Hubble radius, k — aH, at present. 



quantities, as long as they are less than of order unity. In particular, this holds true (to better than 1 per cent) when 
it is assumed that there are no initial scalar field perturbations either in the ZSG, i.e. {S(j)')^ = Sc/)^ = 0, or in the 
UCG, i.e. (Sc/)')'^ = (^^"^ — 0. So the initial condition for the scalar field perturbations does not seem particularly 
important. 

In Fig. 2 we plot the evolution of 5.^ (a) for perturbations with a comoving wavenumber of A: = 0.5 hMpc~^ between 
z = 100 and the present. All the four scalar field models are shown, together with the cases of a critical-density 
Universe and a flat Universe with Og = 0.4. The integration runs were actually started at z — 1100, but all models 
behave very much as the critical-density case up to z — 100. We chose the wavenumber shown as it corresponds to 
a scale which was already well inside the Hubble radius at the beginning of the integrations, so we need not worry 
about specifying the gauge. In fact, comparing the present amplitude of (5^ obtained in the ZSG and in the UCG for 
the various models under consideration, we find that the density perturbations in the two gauges coincide for values 
of k down to about 10"'^ Mpc""'^, and can therefore be regarded as gauge-independent. For smaller k values one 
needs to be careful to specify the gauge used. 

As mentioned in Section III, Fig. 1 represents the distortion in the square root of the matter power spectrum 
suffered since redshift 1100, P{k) oc S^, relative to a flat Universe with fio — 0.4 and no scalar field. Since the full 
spectrum in that model, including all known effects of radiation and neutrinos, is well known, we can therefore use it 
to obtain the expected present shape of P{k) for our four scalar field models. 

We show the power spectra in Fig. 3, where we simply multiplied the expected present shape of P{k) in a flat Universe 
with Qo = 0.4, Qb — 0.016 and either h = 0.55 or h — 0.6, with no scalar field present, by 6'^/6^{Qo — 0.4). This 
procedure is correct since up to z = 1100 the background and perturbation evolution of any of the scalar field models 
is equivalent to that of a cosmological constant model which has the same values for Hq and h as each individual 
scalar field model. We also show in Fig. 3 the expected present shape of P{k) in a critical-density Universe for the 
two values of h considered. We obtained the expected present shape of P{k) for both the critical density and the 
cosmological constant cases by assuming a Harrison-Zel'dovich (scale-invariant) primordial power spectrum for the 
energy density perturbations and using the cold dark matter transfer function initially derived by Bardeen et al. pO[ , 
and later modified by Sugiyama |4l|| to include the contribution of baryons. 

The critical-density and cosmological constant models are normalized to COBE pS]. The amplitude of the other 
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FIG. 2. The evolution in the ZSG of the perturbations with a comoving wavenumber oi k = 0.5 /iMpc~^. The dash-dotted 
and dotted lines represent the cases where respectively the Universe has critical density and the Universe is flat with f2o = 0.4, 
with no scalar field present. Solid lines show the EXP models and dashed ones the PNGB models. 

models, which has been computed relative to the latter, is almost correct. It includes the two main effects — that the 
early time matter power spectrum in low-density models is higher by a factor l/ilo than in critical-density models 
is encoded in the initial conditions through the different initial value of H , and the growth suppression factor from 
the scale factor dynamics is computed in the subsequent evolution. As each model (except the critical-density one) 
has the same redshift of matter-radiation equality, all that is omitted is the line-of-sight contribution to the Sachs- 
Wolfe effect, which requires a full Boltzmann code for accurate computation. The effect of this term on the COBE 
normalization is non-existent for the critical-density case, and known to be negligible in the cosmological constant 
case (see e.g. Rcf. [^^). It may be slightly more significant for the decaying cosmological constant case and 
future observations would probably require an accurate computation, while present ones do not. 

We only show the power spectra for values of k larger than 10"'^ /iMpc~^. On smaller scales the amplitude of the 
power spectrum is gauge-dependent, though of course well-defined in any particular gauge. 



The results presented in this paper extend previous work on the effects of a presently-existing scalar field, with 
regard to the evolution of energy density perturbations in a pressureless fluid. In accordance with general covariance 
we allow for the presence of spatial perturbations in the scalar field. Contrary to previous authors, we explicitly relate 
the energy density perturbations in the pressureless fluid to those in the scalar fleld through the adiabatic condition, 
as expected if those perturbations arose from a period of inflation in the very early Universe. 

The study of the perturbation evolution was performed both in the zero-shear and uniform-curvature gauges, in 
the latter case for the first time. 

We considered two possible potentials for the scalar field, an exponential and that associated with a pseudo-Nambu- 
Goldstone-boson field. Each of these potentials has two degrees of freedom, and we worked with two different sets of 
parameters for each potential. Those models, EXP2 and PNGB2, which were chosen so that the scalar field becomes 
dynamically important only very recently, at z ~ 1, and is presently behaving as a cosmological constant, yield a 
shape for the power spectrum of energy density perturbations in the matter component extremely close to that one 
obtains in a cosmological constant model where fio = 0.4 and h = 0.6, though the amplitude is a few percent smaller. 



V. DISCUSSION 
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FIG. 3. The power spectrum of the energy density perturbations in the matter component for the models, all normalized 
to COBE as described in the text. The line styles represent P{k) for the same situations as in the previous two figures. Note 
that the PNGB2 and EXP2 models give a very similar spectrum to the standard cosmological constant case. 

The situation is rather different for the other two models, EXPl and PNGBl, chosen so that the scalar field becomes 
dynamically important at a much earlier time, around a redshift of 100. They presently behave as pressureless matter. 
In the case of model EXPl the shape of the power spectrum of energy density perturbations in the matter component 
is again very close to that one obtains in a cosmological constant model, this time with fig — 0.4 and h — 0.55. 
Its amplitude is however only about 5 per cent of that for the cosmological constant model. Model PNGBl has 
the further interesting feature that the shape of P{k) is altered for values of k below O.Ol/iMpc"^. Unfortunately 
this corresponds to scales above around 100/i^^Mpc, at the limit of what present cluster and galaxy surveys can 
probe. The feature in the shape of P{k) for model PNGBl is due to oscillations in the rate of decrease of the Hubble 
parameter when the scalar field becomes dynamically important at z 100. On smaller scales the suppression of the 
growth of the energy density perturbations in the matter component is similar to that in model EXPl. 

The results just described suggest that the presence of spatial perturbations in a presently-existing scalar field 
substantially affects the evolution of energy density perturbations in the matter component only if the scalar field 
has contributed significantly to the total energy density in the Universe for several Hubble times. If the scalar field is 
becoming dynamically important only now, the presence of spatial perturbations does not seem to have much effect. 

The similarity between the scalar field models EXP2 and PNGB2, and a flat model with JIq = 0.4, h = 0.6 and no 
scalar field, makes the first two models attractive substitutes of the latter. Further, given that the above mentioned 
flat model is able to reproduce all the most reliable observational data presently available on large-scale structure and 
cosmic microwave background temperature anisotropics we expect models EXP2 and PNGB2 to do as well. As 
they are quite similar to the standard cosmological constant case, a detailed evaluation of the supernova constraint 
would be interesting; Frieman et al. suggest that the constraint will weaken so they should be viable. 

With regard to models EXPl and PNGBl, their much smaller growth for the matter perturbations seems very 
problematic. We expect the dispersion of the density contrast smoothed on spheres of 8 Mpc, usually represented 
by as , required for both models so that they can reproduce the present abundance of high mass X-ray emitting galaxy 
clusters to be the same as for the critical-density case, due to the two models being dynamically equivalent to it 
since well before redshift 10. This conservatively requires ag in the range 0.45 to 0.8 for both models [Q. On the 
other hand, as the growth suppression factor of the matter perturbations is about 4.5 times larger in the EXPl and 
PNGBl models than in a flat model with fig — 0.4 and h = 0.55, to a first approximation this means that the value 
of (Tg implied by COBE for the two scalar field models would be about 4.5 times smaller than for the flat model pl| . 
Inclusion of the integrated Sachs- Wolfe effect may lead to a slight further decrease in erg [|l| . If the primordial power 
spectrum of energy density perturbations is assumed scale- invariant, then this would mean that erg < 0.18 for both 
EXPl and PNGBl One would need a very 'blue' primordial power spectrum, with at least n > 1.45, for cg to 
reach the minimum requirement of 0.45. Ferreira and Joyce advocate a much larger value of flo, which resolves 
the amplitude problem but raises the question of whether the change in the shape of the spectrum can really be 
enough. We are presently carrying out a full investigation of all these models against a range of observational data. 

As in both EXPl and PNGBl the scalar field goes through a period at a redshift of about 100 when it behaves 
dynamically like a cosmological constant, we expect the full angular anisotropy spectrum for the cosmic microwave 
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background radiation in these models to display a distinctive signature. This clearly merits detailed investigation. 



Note: As we were completing this paper, a preprint by Caldwell et al. |23| ] appeared which covers similar issues. 
They use an enhanced Boltzmann code to generate microwave anisotropy power spectra. Their matter power spectra 
appear in good agreement with ours. They did not use adiabatic initial conditions for the scalar field perturbations, 
but we have shown that this is unlikely to have significant effects. 

ACKNOWLEDGMENTS 

P.T.P.V. was supported by PPARC and A.R.L. by the Royal Society. We thank Scott Dodelson, Pedro Ferreira 
and Martin White for discussions, and Wo-Lung Lee for pointing out an error in an earlier version of this paper. We 
acknowledge use of the Starlink computer system at the University of Sussex. 

APPENDIX A: THE EVOLUTION EQUATIONS 
1. Gauge considerations 

A cosmological perturbation is defined by means of a correspondence between an arbitrary background spacetime 
and the real physical Universe. A gauge transformation is a change in this correspondence, keeping the background 
spacetime fixed. Therefore, in general the value of the perturbed part of some physical quantity will not be invariant 
under a gauge transformation. Further, the degrees of freedom due to gauge transformations give rise to spurious 
unphysical modes in the solutions to the evolution of the perturbed part of gauge-dependent quantities, which can 
always be removed by a convenient gauge transformation. In order to avoid these spurious modes one either evolves 
gauge-invariant quantities related to the gauge-dependent quantities one is actually interested in, obtaining the latter 
from the former at any one time, or one has to provide a gauge-fixing condition which completely specifies the way 
through which spacetime is to be split into background and perturbed components. 

The linear analysis of cosmological perturbations was initiated by Lifshitz, with a seminal paper published in 1946 
psf , who used the so-called synchronous gauge-fixing condition. Unfortunately this condition, though considerably 
simplifying the perturbation equations, still leaves a residual gauge degree of freedom. The spurious modes thus 
arising are difficult to distinguish from the real physical ones, and their identification was a source of controversy for 
some time. The use of gauge-invariant quantities in the calculation of the evolution of cosmological perturbations only 
really took off with the paper by Bardeen in 1980 [Q. Though it avoids the problem of obtaining unphysical modes 
in the solutions to the perturbation equations, it really does not offer any extra advantage over the gauge-specific 
methods which remove any gauge degree of freedom by completely fixing the background/perturbed splitting. 

In this paper we will use gauge-specific methods to solve the perturbation equations. We will adopt the notation 
and use the equations laid down in a series of papers by Hwang 1 29 - 33j . We will be solely interested in the evolution 
of density (scalar) perturbations. The system is composed of an ideal fluid plus a single minimally coupled real scalar 
field, 0, evolving in a background Einstein-de Sitter Universe. Given that the spatial part of the background spacetime 
is thus homogeneous and isotropic, the perturbations in any physical quantities will necessarily be gauge-invariant 
under purely spatial gauge transformations. We will therefore only worry about the temporal gauge transformation. 

2. Notation and general equations 

We will now introduce the gauge non-specific perturbation equations obtained by Hwang. They relate the perturbed 
part of the metric variables, a (perturbed part of the lapse function), Lp (perturbed part of the spatial curvature), x 
(perturbed part of the shear) and k (perturbed part of the expansion scalar), to the perturbed part of the matter 
variables, e = e-y-l-e,/, (perturbed part of the total energy density), vo = vo-si + vo^ (perturbed part of the total pressure)^ 
and = -|- (perturbed part of the total energy density flux, or total fluid four-velocity, depending on the frame 
chosen). We have 



^Here we changed the notation from -n to -co to avoid confusion with the number tt. 
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^7 = (7 - l)e7 ' (Al) 

and one can derive that 

e^^<j>d<p-<j>'^a + V,^S(t>, (A2) 

^ 4'S(j) - (pa -V,^S(I) , (A3) 
= (A4) 

where 6(j) is the perturbed part of the scalar field. We will be particularly interested in the evolution of the fractional 
perturbation in the energy density of the ideal fluid, S-y = e-y//i-y. 

We will express the perturbed parts of both the metric and matter variables by means of Fourier expansions. For 
example, 

5</)(x,t)=^<50kWe*-% (A5) 

k 

where 

SMt) - (5<^(x,t)e-* '^dx, (A6) 

being /c = |k| a fixed comoving wavenumber. The Fourier expansions are made in a large enough box that the induced 
periodicity is irrelevant. 

As in this paper we are only interested in the linear evolution of cosmological perturbations, we will assume that the 
different Fourier modes for each variable behave independently of each other. We will drop the suffices k identifying 
each Fourier mode in order to lighten the notation. 

In the case of our system, formed by an ideal fluid plus a minimally coupled real scalar field, the perturbation 
equations take the form [ p9| , p0p^ ], 

^2 

Sip = SHa -K+—X, (A7) 

-—ip + HK^ -47rG(e^ + (l>S(t> - (j)^a + ¥,4, Scj)) , (A8) 

K-^X = -127rG(*^-(/)(5(/)), (A9) 
± + Hx = a + ^, (AlO) 

k + 2Hk = f ^ - 3h] a + 47rG[(37 - 2)e^ + 4050 - A^a - 2V,^ 5(1)] , (All) 



fc2 

fi-y + iH^e-y = ^^i-y{K-'iHa) + —-^-f, (A12) 

+ SiJI*^ = - (7 - l)e^ , (A13) 

54> + mScf, + {^ + y,4>4> 0(k + a) - (3i?0 + 2V,^ )a . (A14) 

It should again be stressed that these equations were obtained without reference to any gauge-fixing condition. No 
anisotropic pressure term appears as in both the case of an ideal fluid and a minimally coupled real scalar field the 
anisotropic pressure is zero. The last three equations are, in order, the energy and momentum conservation equations 
for the perturbations in the ideal fluid, and the energy conservation equation for the perturbations in the scalar field. 
The momentum conservation equation for the perturbations in the scalar field is identically satisfied. 

The most obvious and fundamental gauge-fixing conditions follow from requiring that the perturbed part of one of 
the metric or matter variables is zero. We thus have: the synchronous gauge, a = 0; the uniform-curvature gauge, 
(/? = 0; the zero-shear gauge, % = 0; the uniform-expansion gauge, k = 0; the uniform-density gauge, e = 0; the 
uniform-pressure gauge, 137 = 0; and the comoving gauge, VE" = 0. Except for the synchronous gauge, all the other 
gauge-fixing conditions completely remove the gauge modes from the solutions to the perturbation equations. We will 
use two of these gauge-fixing conditions to derive two different sets of perturbation equations from the system given 
above. With the aid of expressions relating quantities in the two gauges we will thus be able to estimate the numerical 
errors arising from the numerical integration of both sets of perturbation equations. We will consider the zero-shear 
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gauge (ZSG, also known as the longitudinal or conformal Newtonian gauge [0), and the uniform-curvature gauge 
(UCG). These choices are the ones which lead to the two simplest sets of perturbation equations, thus decreasing the 
probability of numerical errors creeping into the solutions. The two sets can b e obtained by simply getting rid of x 



and X in the case of the ZSG, and ip and ip for the UCG, in Eqs. (A?) to ( A14 ). 



APPENDIX B: NUMERICAL SOLUTIONS 

1. Background and perturbation equations 

The evolution of the background variables H and (p is obtained by numerically solving the system of first-order 
differential equations formed by Eq. (||) and the two first-order differential equations that can be obtained from Eq. (Q), 

^^ = -4,G(ai//V^,„). (B3) 

where /i^ is given by analytically solving Eq. (j^), 

The suffix '0' indicates present-day values as usual. Note that the independent variable has been changed from 
coordinate time, to the scale factor, a. They are related by the first-order differential equation 

# = a-^H-^ , (B5) 
da 

the integration of which gives the time elapsed in the Universe between two given values of the scale factor. Derivatives 
with respect to the scale factor will be represented by a prime. We choose the independent variable to be the scale 
factor as it is easier to work with numerically and is more meaningful from the point of view of structure formation. 
We have 8 perturbation equations for 11 dependent perturbation variables in both the ZSG and the UC G: e^, e' 
Scj), Scj)' ^ 5(f)" ^ a and n in either, along with ip and (p' in the ZSG, or x and x' in the UCG. Eqs. (A7) to ( A1C| ), 



( A12| ) and ( A13| ) will be used to describe the evolution of a, k and the quantities associated with the ideal fluid, e^. 



'^-y and W^iu terms of 50, 54>' , Lp and Lp' (ZSG) (the last two variables are replaced by x and x' in the UCG), 
and the background variables. We thus have in the ZSG, 

a = -p,, (B6) 
-ZHlp - 3aHp' , (B7) 

p P 3H^(p + 3aH^(p' 2„2<^XA' 2rr2,/2 r, mo^ 

H ——^ aH(pd<p-aH(j} ip -V^^ocp, (B8) 



K 



^1 



4ttG a2 AttG 



*. = ^^i^+«W0, (B9) 



AttG 
e 



-37^-37/i^(^'-t-3^, (BIO) 



a2 aH 



and for the UCG, 



a = Hx + aHx' , (B12) 



K = ( 3i/^ + ^ ) X + 3aH\' , (B13) 
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^7 


/ T O \ , , — ^ -_--_-Q / 

V a^y ^-kG 47rG 




*7 = 








-37— + 7Ai7 — ^ + — ^ ' 
a ail a a/i 






7^7 (7 1) jj- ■ 

a a a_H 





(B14) 

(B15) 
(B16) 
(B17) 



Using these expressions we can now convert Eqs. (All) and (A14) into the following second-order differential equations, 



H' 



a H 

(37 + 2) 



= - - + ^ - ^ + K 



fc2 



IT 
ll 



37 

n2 



9H' P 

-FF + (7-l)^-4^G(7-2) 

atl a^H'^ 



AttG 



(7 - 2)(l)'d(j)' + 7 



(B18) 



(B19) 



in the ZSG and 



- + ^ + 4^G(7 - 2)a0'2 
a 



fc2 
«2 



+47rG7a</)'y,, 



2H'(h' 



aH 



- 47rG(7 - 2)a^H(j)' 



/3 



a a^ii a^ ail 



X. 



x" = 



37 



— + 3^ - 4^G(7 - 2)a0'2 
a H 



^ + S + (7-l)^-47rG(7-2)0'2' 
aii a^ a^H^ 



-47rG 



(7 



H 



+ 7- 



l2if3 



(B20) 



(B21) 



in the UCG. Each of these equations can be split into two first-order differential equations, in the same way as we did 
for Eq. (jj), which we will then numerically integrate in order to determine the evolution of 50, Sc/)' , Lp and ip' (ZSG) 
(X and x' in the UCG). 

In all we will need to simultaneously n ume rical ly in tegrate a system of seven first-order ordinary differential equa- 
tions, formed by the backgroun d Eq s. ( Bl) to ( |B3| ), where expre ssion (B4) gives fij, and the four perturbation 
equations resulting from either ( pl8| ) and (|b19D (in the ZSG), or ( p20| ) and ( p2l| ) (in the UCG). The question of 
initial conditions for this procedure is addressed in the main text of this paper. 



2. Relations between quantities in different gauges 

Once we have calculated the evolution of the perturbation variables in s ome part icular gauge, we can use gauge- 



invariant variables and the gauge non-specific set of perturbation Eqs. (A7) to (A14) to obtain the evolution of such 
variables in any other gauge. 

We will use this possibility to control the errors arising from the numerical integration of the perturbation equations. 
We will express both our initial conditions and the final results for the perturbation variables in the ZSG, and use the 
UCG simply as an estimator of the numerical integration errors. As an example, we will derive the relations between 
the perturbations in the energy density of the ideal fluid in the ZSG and the UCG. The change from the ZSG to the 
UCG, and vice-versa, for the other perturbation variables can be obtained in a similar way. 

The quantities 

e^ = e^+3H{fi^+p^)x, (B22) 
^^ = -^^ + i^^+p^)x, (B23) 

and 
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H 



(B24) 
(B25) 



are invariant under temporal gauge transformations, as can be seen by using the relations provided by Hwang [ pO[ . 
The first two variables simply take the values of e-y and ^'-y, respectively, when these quantities are calculated in the 
ZSG, while the same occurs for the last two variables with relation to the UCG. We thus want to express in the UCG 
as a function of and VE"^, in order to obtain the initial value of to be used in the UCG calculations from that 
originally given in the ZSG. Also, we want to know how to obtain from e"^ and so that we can compare the 
final value of obtained in the two gauges. 
In the UCG we have 



^7 ^7 



and, by using Eqs. ( [AS| ) and ([A9|), 

fc2 



Through some algebraic manipulation of the above relations we then obtain 

„2 
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where and 4'^ are defined in an analogous fashion to the ideal fluid gauge-invariant variables. 
On the other hand, in the ZSG we have 

= e!^ - 3{^iy+p^)if, 
and, by using the same equations as previously. 



Hence, we get 



(B26) 



(B27) 



(B28) 



(B29) 



(B30) 



(B31) 



where again and ^'^ are defined in the same way as the ideal fluid gauge-invariant variables. 
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